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A scheme of left-handed metamaterial (LHM) composed of superconduct- 
ing quantum interference devices (SQUIDs) and conducting wires is proposed. 
The permeability of a probe field can be smoothly tuned over a wide range 
with another electromagnetic (coupling) field due to quantum interference ef- 
fect. Similar to electromagnetically induced transparency (EIT) of atomic 
systems, the absorption of the probe field can be strongly suppressed even 
in the case of negative permeability. There are two passbands of negative 
refractive index with low loss, which can be tuned with the coupling field. 

PACS numbers: 42.50.Gy, 78.20.Ci, 03.67.Lx, 42.50.Ct, 32.80.Qk 

Recently, research on left-handed materials (LHMs)pQ has attracted con- 
siderable attentions. The LHMs are a kind of metameterial which have nega- 
tive permittivity and negative permeability , which lead to negative refractive 
index in a narrow frequency band. Many novel phenomena can occur on LHM 
such as superprism, perfect flat lens, inverse light pressure, and reverse Doppler 
and Vavilov-Cherenkov effects . Large nonlinearity can also occur, e.g. the 
bistable transition of permeability from positive to negative 2 j. The quan- 
tum phenomena in LHMs have also attracted attentions such as the modified 
spontaneous emission of atoms in LHM [31 II]- However, the background of the 

ordinary LHM such as array of conductor lines and slit-ring resonators (SRRs) 
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is a 'classical' system where the negative e and negative /i arise from the clas- 
sical plasma oscillations. In this sense, the ordinary LHMs are quite different 
from quantum systems such as atomic gases. On the other hand, in ordinary 
situations, the magnetic response of atomic gases to a laser field is too weak 
to generate negative permeability. How to enhance the magnetic response of a 
quantum system is essential to realize LHMs at optical frequency band. There 
may be other ways to obtain negative permeability, e.g., recently, a scheme 
of electromagnetically- induced left-handedness in atomic gases is proposed 
where large density of atoms is necessary and special restrictions on frequencies 
of driven fields and atomic transitions are required. In contrast, an artificial 
quantum system with high magnetic sensibility can has much more advantage 
than atomic systems to realize negative permeability. 

In this letter we propose a new kind of LHM composed of superconducting 
rings with Josephson junctions and conductor wires. For simplicity, the wires 
are assumed to be normal conductors. It should be note that quite different 
from our work, recently analysis and design of superconducting transmission 
lines are presented by Salehi et al [H], and more recently, Ricci et al have ex- 
perimentally researched on a metamaterial that employs superconducting Nb 
metals and low-loss dielectric materials, in which case negative effective index 
passband are seen between 50MHz to 18GHz [7j. The focus of our work is on 
the negative permeability which arises from quantum feature of the compos- 
ite. In contrast with the 'classical' SRR, in our model the split of the ring is 
replaced by a Josephson junction, which is essential to the quantum feature 
of the LHM. We will show analytically that due to the quantum interference 
effects |S] the permeability can be tuned over a wide range with an external 
microwave field. Also, we will show that the absorption of the medium for the 
probe field can be strongly suppressed even in the case of negative refractive 
index. 
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A scheme for supper conducting LHM is show in Fig.l, where the compos- 
ite is composed of superconducting rings with Josephson junctions, and, an 
effective electric medium with effective electric permittivity e as a background, 
which, e.g., can be an array of normal conducting wires. The radius of the 
ring is denoted by a , and the period of the array is denotd by d. Schematic of 
the potential energy and the first six eigen energies of the SQUID is shown in 
Fig. 2. We assume a probe microwave field is interacting with the composite. 
The wavelength of the field is assumed to be much longer than the period, 
i.e., d « A. For a classical cylindrical SRR system, the permeability \i can 
be given by the relation (according to Ref.[2]) 

B(w) = H a {u) + FH'{u), (1) 

where H x (u>) is the alternating external magnetic field and H'{uj) is the addi- 
tional magnetic field induced by H x (tu), which determines the magnetization 
of the composite, and F = -Ka?/d 2 is the fraction of the structure with a 
being the radius of the ring and d being the periodicity of the array. Also, 
a « d « A is assumed. 

The permeability at angular frequency u> can be given by 

„ M = 1 + f|£> (2) 

where 4>{oj) is the flux of frequency uj induced by the external microwave field, 
i.e. 4>{oj) = H's with s = ira 2 being the area of the ring. For the quantum 
LHM here, the SRRs are replaced by the SQUIDs, then (ft(ui) becomes an 
operator although the external driven fields are assumed to be classical. In 
order to calculate fj,, cp(uj) in Eq. © should be replaced by the quantum 
averaging of it, i.e. < 4>(uj) >■ 

The Hamiltonian of a SQUID (a ring with a Josephson junction) can be 
given by[S] 
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with the potential of the SQUID being 



V(x) = -muj 2 LC (x - x') - —^mu} 2 LC [3cos{2Ttx) (4) 

where x = <fi/(f>o, m = C(f)Q, lo 2 lc = [3 = 2ttLI c /4>o, and x' = 4> x /4>q. 
Here 4> is the total flux in the ring, L is the ring inductance, 4> x is an external 
applied magnetic flux to the SQUID, I c is the critical current of the junction, 
C is the capacitance of the junction, and (j>o(= h/2e) is the flux quantum. 
We consider a realistic SQUID system which can be described by use of the 
parameters as in the work of Zhou etal [Sj, where L = WOpH, C = 40/.F, and 
I c = 3.95fiA, leading to lolc = 5 x 10 n rad/s, and f3 = 1.2. The external DC 
magnetic field parameter x' is taken to be —0.501. 

The interaction between the SQUID and microwave fields, which are as- 
sumed to be linearly polarized with their magnetic field perpendicular to the 
plane of the SQUID ring, is described by the time dependent potential 

Vint{x, t) = m ^Lc( x ~ x'){ecosu)t + £ c COSU) c t). (5) 

In the interaction picture ,the dynamics of the system is governed by the 
Schrodinger equation 

ifcj^(t) >=V int \m > . (6) 

Here e, e c are microwave magnetic flux of the probe field and the coupling 
field in units of cpQ. The frequency of the probe field and the coupling field are 
chosen near resonant with transition |0 > — \n > and |1 > — \n >) respectively, 
where \n > denotes the n-th eigen state of Hq. This is usually referred to as 
the three-level approximation. In this case, the wave function can be written 
as 

\ip(t) >= co(i)|0 > +ci(t)|l > +c n {t)\n >, (7) 
and the interaction Hamiltonian after the rotating-wave approximation, can 
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be written as [S] 

(8) 



Vim = H^e iSt \0 >< n\ + n c e- iSct )\l >< n\ + h.c.-i^-\l >< l|-iy|2 >< 2|). 



Here the Rabi frequency f2 and f2 c are denned as 

^ = -xo n mojl c e/h, fi c = -xi n mu? LC e c /fi, (9) 

where xq u and xi n are defined as < 0\x\n > and < \\x\n > respectively. Also, 
here 5 = to — u>o n , 5 C = to c — u; ln , 71 and 7 n are the background decay rate of 
the state |1 > and \n >. The equations for amplitudes Co, c\, and c n can be 
easily obtained from Eq.® as 

■dco _ n„ 
1 dt ~ 2 Cn 

i<kg. = _ (5 + i^ )Cn + q Co + & C1 , (10) 



dt — V 2 / ^ 2 2 

^ = -(' 5 -^ + i T) c l + ^ c - 
We assume that 71 << j n , where 7„ can be taken to be about 1GHz according 
to Ref. |Hj . We assume that initially the system is prepared in the ground 
state |0 >, and for simplicity, we assume the probe field is weak enough, i.e. 
^ << 7n 5 ^ << ^c- In this case the steady-state solution can be obtained as 
in ordinary EIT systems. 

It is easy to prove that the mean value of the flux (< c/)(uj) >) can be given 

by 

<m>—.&., (id 

where a = — L £ °" = j^Xq u . Here a is a parameter measures the sensibility of 
the effective medium to the probe field, which has the dimensions of frequency 
(Hz). For n = 4, x 4 = 5.39798 x 1CT 3 , which is given by Ref . [Ill El • Taking 
L = WOpH, we obtain a = 11.815GHz. This value will be used in the following 
analysis. 

We first consider the case of O c = 0, i.e. of two-level system (non-EIT 
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system). According to Ea.l|2T). QllJI and f|X0|> . the permeability can be give by 



, It is clearly a Lorentz profile. When the probe field is near resonant with the 
transition |0 > —\n >, strong absorption will occur. The absorption becomes 
sharp with the decrease of 74. Similar to classical SRR composite, negative [x 
can occur (see Fig. 3(a)) if a is large enough. The condition for the negative 
H can be easily deduced from Eq. ^]as that Fa > 274. Also, it is clear that 
there is only one frequency band of negative Re(/x). However , it is worth to 
note that the //-spectrum for the quantum composite here is similar but not 
the same as that in classical composite 2 . 

If the coupling field is applied to the system (Q c ^ 0), the magnetic 
response could be significantly modified. For a strong coupling field, i.e. 
f2 c >> f2, 74, from Eg. ((2)) . (| 11 1) and (|1U|). the permeability can be given by 



The spectrum of \i in three-level case is show in Fig. 3(b). It can be easily 
seen that the transparency occurs when 5 = 5 C = 0, which is similar to EIT of 
atomic system except that here the coupling field is used to control the per- 
meability fj, instead of the permittivity e of the probe field. Superconductive 
analog to electromagnetically induced transparency that utilizes superconduc- 
tive quantum circuit designs of present day experimental consideration has 
also been investigated by Muralial et al [S]. The width of the transparency 
window can be defined by the frequency distance between the two absorption 
peaks which is about f2 c due to Autler-Townes splitting. The frequency band 
for low absorption and large refraction can occur over a wide range for large 
Q c . On the contrary, the group velocity of the probe field can be reduced to 
near zero as in atomic EIT systems. If Fa is very small, /x is positive for all 
probe detunings (<5), which can be easily seen from Eq. (|13|) . If Fa is large 



F a 
~2 5 + ij 4 /2 



(12) 



fi = l- 




(13) 
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enough, however, negative Re(/u) can occur near the transparency window. 
The condition for negative /i can be deduced from Eq. as that 

Fa > 2 74 . (14) 

This condition is the same as that for nonEIT case, but the frequency band 
for negative \x in the EIT case and nonEIT case is quite different. In the EIT 
case there are two frequency bands for negative fx, which can be given by 

f-ix-) < S / 7 4 < f-(x + ); f + {x-) < 5/74 < /+(z+), (15) 

where 



and 



f±(x) = V - , (16) 



x± = -{g±yj<?-l) (17) 

with 5 = g(> 1) 

For example, we consider the case of a = 11.815GHz as above. If taking 
F = it x 0.14, we can obtain that Fa/2 = 2. 5982GHz. In this case negative 
Re(/x) can occur if the condition 74 < 2. 5982GHz is met. 

It should be note that, when 5 = 5 C and 71 = 0, exactly zero absorption 
occcur (Im(/i) = 0), but in this case Re(/i) = +1. When ReQu) is negative, 
Im(/i) is always nonzero, i.e. negative refractive index without loss can never 
occur. However, the absorption can be strongly suppressed near the condition 
of 8 = 5 C . It is similar to ordinary EIT in atomic systems, but here the 
zero absorption occurs due to magnetic response instead of electric response, 
and, here, the magnetic is so strong that negative Re(/i) can occur near EIT 
condition. The permeability (/i) can be smoothly tuned over a wide range with 
the coupling detuning 5 C (see Fig. 4) and its Rabi frequency O c (see Fig. (5)). It 
is worth to note that the tunability is due to the quantum interference effect 
which is absent in classical normal conductor systems. The permeability /j, 
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and refractive index n can be smoothly tuned by the coupling field over a 
wide range, and the continuous transition from Re(/x) = 1 to Re(^t) = — 1 is 
possible. Prom the /x-spectrum it is easy to see that there are four frequencies 
for Re(//) = — 1, but only two of them are for low loss (small Im(^)), one of 
which is in EIT window, while the other one is out of EIT window. It is should 
be note that both of them can be tuned with the coulpling field. For this two 
frequencies, Im(/i) — > when 74 — > and 71 — ► 0. 

In order to calculate the refractive index, the permittivity e and perme- 
ability fM can be written as 

e= | e |e^ e ,/i= |/i|e^ m , (18) 

where < 4> e < it, < 4> m < it. Then the refractive index n can be given by 
3 

n = y^He^^™)/ 2 (19) 

For the sake of simplicity, we consider the case where the effective electric 
permittivity e is generated by an array of normal conducting classical current 
wiresjHlH]. In this case 

e = l+ a I 6 . (20) 

where Lo pe and wy e are the parameters which measure the effective plasma os- 
cillations and 7 e measures the loss of the wires. The refractive index spectrum 
of the composite is shown in Fig. (6), where (a) is for e(<5), (b) is for /J,(S), and 
(c) is for n(5). For simplicity, we consider the frequency band of negative e 
which is far from resonance, in which case e is slowly varing with 5, and Im(e) 
is small. We find in this case there are two minima on the Im(n)-spectrum. 
The absorption minima arise from two reasons: (1) small Im(e) and small 
Im(/i); (2)Re(/i) and Re(e) are simultaneously negative. On the contrary, if 
Re(/x)i2e(e) < 0, Im(n) will be large even if both of Im(ju) and Im(e) are in- 
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finitely small. Two passbands of negative refractive index appear around the 
two minima, both of them are sensitive to the coupling field. 

In conclusion, we have analytically investigated the metamaterial com- 
posed of superconducting rings with Josephson junctions and conducting wires. 
It is found that negative permeability for a probe microwave field can occur 
when g = ^ > 1 . There are two passbands of negative refractive index with 
low loss, which can be tuned with the coupling field. To our knowledge, it 
is for the first time that the quantum left-handed metamaterial being com- 
posed of SQUIDs has been investigated. One of the advantages of this kind 
of LHM is that the negative refractive index with low loss is easy to obtain 
due to that the permeability can be smoothly tuned over a wide range includ- 
ing fi = ±1 + ie (e — > 0+). In the quantum composite, some new physics 
could be found, such as that associated with transient properties of the left- 
handedness, large nonlinearity due to the quantum interference effect and that 
due to strong magnetic response of the composite. Also, the tunability of the 
composite via quantum interference effect could be used in quantum informa- 
tion process |13j . On the other hand, the LHM with SQUIDs can be regarded 
as a bridge between the 'classical' LHMs (e.g. composed of SRRs and wires) 
and the 'quantum' systems composed of microscopic particles(e.g. atoms or 
quantum dots). Hence, further research on the quantum LHM could facilitate 
the realizing of the tunable LHMs of high-frequency (towarding to optical) 
band. 
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Figure Captions 

Fig. 1. Schematic of the composite metamaterial structure composed of 
superconducting rings with Josephson junctions (SQUIDs) and the effect 
electric medium such as conductor wires. Each SQUID is an artificial 
A-configuration three-level system which coupled with a probe microwave 
field and a coupling microwave field. 

Fig. 2. Schematic of the potential energy and the first six eigen energies of 
the SQUID. The energies of the ground state |0 >, meta-stable state |1 >, 
and the excited state |4 > are 7.81984mev, 7.90183mev, and 8.14057mev, 
respectively. 

Fig. 3. Real part (solid curve) and imaginary part (dashed curve) of the 
permeability (/z) versus the probe detuning S for 5 C = 0, where F = ir x 0.14, 
74 = 0.0423, 7i = O.I74, (3 = 1.2, x' = -0.501. (a) is for n c = and (b) is for 
S7 C = 0.5078. All parameters are in units of a, which is taken to be 
11.815GHz. 

Fig. 4. Real part (solid curve) and imaginary part (dashed curve) of the 
permeability ([/,) versus the coupling detuning 8 C for the case of 8 = 0. The 
parameters are the same as that in Fig. 3 

Fig. 5. The band edges of negative Re(/z) versus the Rabi frequency of the 
coupling field Q c . Negative Re(/i) occurs when S is between the two solid 
curves (bandl) or that between the dashed curves(band2). The parameters 
are the same as that in Fig.3. 

Fig.6. The spectrum of permittivity (a), permeability (b), and refractive 
index(c), where u> pe = 1.5cjo4, wtc = 0.43wo4, 7 e = 0.1wo4- The solid line is 
for the real part of them, while the dashed line is for the imaginary part of 
them. Here W04 = 48.727GHz. Other parameters are the same as that in 
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Fig.3. 
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